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1  Introduction 


These  notes  contain  the  slides  of  the  course  Unifying  Parametric  and  Implicit  Surface  Repre¬ 
sentations  for  Computer  Graphics,  given  at  SIGGRAPH  90.  They  begin  with  commentaries 
and  references  that  introduce  the  material  to  individuals  who  have  not  attended  this  course. 

Of  necessity,  the  subject  is  mathematical,  using  many  concepts  from  elementary  algebraic 
geometry.  The  concepts  do  have  an  intuitive  interpretation  that  can  guide  the  reader  using 
those  concepts  even  though  he  or  she  may  not.  be  aware  of  some  of  the  finer  points.  Those 
intuitions  are  only  sketched  here,  but  can  be  found  in  much  detail,  for  instance,  in  [13],  in 
Chapters  5  through  7. 

The  material  is  organized  as  follows:  After  reviewing  some  basic  facts  from  mathematics, 
methods  for  parameterizing  implicit  curves  and  surfaces  are  presented,  concentrating  on  how 
to  deal  with  monoids.  Monoids  allow  an  especially  easy  conceptual  approach,  but  one  can 
appreciate  some  of  the  technical  complexities  when  studying  the  method  in  the  case  of  cubic 
curves. 

Thereafter,  we  discuss  how  to  convert  from  parametric  to  implicit  form.  Many  authors 
have  discussed  these  techniques,  and  only  a  very  limited  perspective  is  developed  here.  A 
broader  description  of  the  subject  is  found  in  [13].  in  Chapters  5  and  7.  and  also  in  [15], 
where  the  problem  of  faithfulness  is  conceptualized  and  discussed.  In  particular,  resultants 
are  covered  in  detail  in  the  theses  by  Sederberg  [17]  and  by  Otyionh  [9].  We  do  not  discuss 
multivariate  resultant  formulations,  but  refer  to  [9]  for  further  reading  on  the  subject. 

Conversion  between  implicit  and  parametric  form  is.  in  general,  an  expensive  computa¬ 
tion.  It  is  therefore  worth  considering  alternatives.  One  such  approach  is  to  view  parametric 
curves  and  surfaces  as  manifolds  in  higher-dimensional  spaces.  Such  a  view  no  longer  has 
to  distinguish  between  implicit  and  parametric  representations,  and  the  methods  it  develops 
apply  to  both  equally  well.  We  will  not  discuss  this  approach  here,  and  the  reader  is  referred 
for  details  to  [1*1.  15]. 

Basic  Mathematical  Facts 

Recall  that  a  plane  parametric  curve  is  defined  by  two  functions 

=  ln(s) 

y  =  M-s) 

and  that  a  parametric  surface  is  defined  by  three  functions 

.v  =  li\(s.t) 
y  =  h'AsG) 

~  =  /'a! -M) 
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We  can  think  of  a  parametric  curve  as  a  map  from  a  straight  line  with  points  s  to  a  curve 
in  the  (x,  y)-plane.  and  of  a  parametric  surface  as  a  map  from  a  plane  with  points  (s,t)  to 
a  surface  in  (x.i/,r)-space. 

The  functions  /i*  will  be  polynomials  or  ratios  of  polynomials  in  s  and  t.  Accordingly, 
we  speak  of  integral  or  rational  parametric  curves  and  surfaces  whenever  the  distinction  is 
critical.  Ordinarily,  the  curves  or  surfaces  are  restricted  in  the  literature  to  a  domain ;  e.g., 
to  the  interval  [0,  lj  or  to  the  square  (0. 1]  x  [0, 1],  Here,  we  do  not  so  restrict  parametric 
curves  and  surfaces. 

Typically,  the  functions  lik  are  presented  in  a  particular  basis;  for  instance,  in  the 
Bernstein- Bezier  basis,  and  this  allows  relating  the  coefficients  of  the  functions  lik  with  an 
intuitive  understanding  of  the  shape  of  the  curve  they  define.  A  suitable  basis  also  affords 
a  wealth  of  techniques  for  combining  patches  of  parametric  curves  or  surfaces  into  larger 
surfaces,  and  to  modify  the  shape  of  the  larger  surface  locally  or  globally,  in  an  intuitive 
manner.  See.  e.g.,  [5,  11]. 

An  implicit  curve  is  defined  by  a  single  equation 

/(•>'».V)  =  0 

and  an  implicit  surface  is  defined  by  a  single  equation 

f(x.y,z)  =  0 

Thus,  the  curve  or  surface  points  are  those  that  satisfy  the  implicit  equation,  so  that  we  no 
longer  think  of  curves  and  surfaces  as  the  result  of  a  mapping.  We  will  restrict  the  function 
/  to  polynomials. 

Since  we  restrict  the  Ilk  to  polynomials,  or  ratios  of  polynomials,  and  restrict  the  / 
to  polynomials,  we  are  dealing  with  algebraic  curves  and  surfaces.  Algebraic  curves  and 
surfaces  include  virtually  all  surfaces  studied  and  used  in  geometric  and  solid  modeling,  and 
in  computer-aided  geometric  design.  Algebraic  geometry  provides  us  with  the  following  key 
facts  about  algebraic  curves;  e.g.,  [20]: 

Every  plane  parametric  curve  can  be  expressed  as  an  implicit  curve.  Some,  but 
not  all.  implicit  curves  can  be  expressed  as  parametric  curves. 

Similarly,  we  can  state  of  algebraic  surfaces 

Every  plane  parametric  surface  can  be  expressed  as  an  implicit  surface.  Some, 
but  not  all,  implicit  surfaces  can  be  expressed  as  parametric  surfaces. 

This  means,  that  the  class  of  parametric  algebraic  curves  and  surfaces  is  smaller  than  the 
class  of  implicit  algebraic  curves  and  surfaces.  There  is  even  a  rigorous  characterization 
of  what  distinguishes  a  parameterizable  algebraic  curve  or  surface  from  one  that  is  not 
parameterizable.  Roughly  speaking,  a.  curve  is  parameterizable  if  it  has  many  singular 
points:  that  is,  many  points  at  which  the  curve  intersects  itself  or  has  cusps.  We  will  not  go 


into  those  details,  because  the  characterization  is  very  technical,  and  the  computations  that 
would  be  needed  to  test  whether  a  curve  or  surface  is  parameterizable  are  quite  complex  and 
time-consuming;  [3], 

We  will  discuss  special  curves  and  surfaces  that  can  be  parameterized  fairly  easily.  They 
include  the  following  cases. 

•  All  conic  sections  and  all  quadratic  surfaces  are  parameterizable. 

•  Cubic  curves  that  have  a  singular  point  are  parameterizable. 

•  Monoids  are  parameterizable.1 

All  parametric  curves  and  surfaces  have  an  implicit  form,  and  we  will  discuss  several 
approaches  for  finding  the  implicit  forms.  However,  it  is  possible  that  a  parametric  surface 
does  not  contain  certain  points  found  on  the  corresponding  implicit  surface.  Some  of  the 
missing  points  can  be  recovered  by  considering  the  surface  parameterization  projectively,  but 
not  all  missing  points  can  be  so  recovered.  Except  for  certain  special  cases,  the  conversion 
between  implicit  and  parametric  form  is  expensive,  and  one  does  not  invoke  the  conversion 
algorithms  lightly. 

Ordinarily,  we  deal  with  affine  spaces  in  which  points  may  be  fixed  using  Cartesian 
coordinates.  On  the  affine  line,  a  point  has  the  coordinate  (xi);  on  the  affine  plane,  a  point 
has  the  coordinates  (xi,X2);  and  in  affine  space,  a  point  has  the  coordinates  (xi,X2,X3). 
contrast,  projective  spaces  add  another  coordinate  x0,  and  consider  a  point  defined  by  the 
ratio  of  its  coordinates.  For  example,  the  projective  line  has  points  (xo,xi),  and  for  all  A  ^  0 
both  (xo,xi)  and  (Axo,Axi)  are  the  same  point.  The  coordinate  tuple  (0,0)  is  not  allowed. 

The  projective  line  can  be  visualized  as  the  pencil  of  lines  through  the  origin,  embedded 
in  the  affine  plane,  as  shown  in  Figure  1.  Here,  the  projective  point  (s,t)  corresponds  to  the 

1  Monoids  are  defined  later. 
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Figure  2:  Parameterizing  the  Unit  Circle 


line 

x0  =  .sA 

.Cl  =  t\ 

To  each  affine  line  point  (xi),  there  corresponds  the  projective  line  point  (l,xi).  Visually, 
the  affine  line  is  viewed  as  the  points  obtained  by  intersecting  the  pencil  with  the  line  x0  =  1. 
The  projective  line  lias  one  additional  point,  with  coordinates  (0, 1),  corresponding  to  the 
line  xo  =  0.  This  point  is  said  to  be  at  infinity. 

The  projective  plane  has  the  points  (xo,xi,X2),  (0*0*0)  is  not  allowed.  For  all  A  ^  0, 
(xo,xi,X2)  and  (Ax0,  Axj,  AX2)  are  the  same  point.  Again,  one  may  embed  the  projective 
plane  into  affine  3-space  by  considering  the  projective  plane  as  the  bundle  of  all  lines,  in 
3-space,  through  the  origin.  The  affine  plane  is  a  subset,  obtained  by  intersecting  the  bundle 
with  the  plane  x0  =  0.  The  additional  points  correspond  to  the  pencil  of  lines  through  the 
origin  that  lie  in  the  plane  xo  =  0,  and  form  the  line  at  infinity.  The  coordinates  of  these 
points  are  of  the  form  (0,Xi,X2). 

Projective  geometry  leads  to  simplifying  many  theorems  by  eliminating  special  cases. 
For  our  purposes,  projective  curve  parameterization  has  the  advantage  that  all  points  on  a 
parametric  curve  can  be  reached  with  finite  parameter  values  without  exception.  This  will 
not  necessarily  be  the  case  for  projectively  parameterized  surfaces. 

2  Parameterization 

2.1  Plane  Algebraic  Curves 
A  Geometric  View  of  Parameterization 

The  geometric  idea  underlying  parameterizing  a  plane  algebraic  curve  is  illustrated  by  the 
unit  circle,  x2  +  y1  -  L  =  0.  We  pick  one  point  on  the  circle,  say  pi  =  (-1,0),  and  consider 


all  lines  through  this  point.  A  line  through  pi  lias  the  equation  y  =  tx  +  (,  where  t  is  the 
intercept  with  the  y-axis.  The  lines  are  indexed  bv  t;  that  is,  to  each  value  of  t  corresponds 
a  specific  line  through  p\.  Each  line  intersects  the  circle  in  pi  and  in  one  additional  point 
p(t).2  The  coordinates  of  p(t)  are  obtained  in  three  steps: 

1.  Substitute  tx  +  t  for  y  in  the  circle's  equation,  obtaining  the  equation 

x2(  1  +  /2)  +  2t.2x  +  t2  -  1  =  0 

2.  Solve  this  equation  for  x.  obtaining 


Observe  that  aq  is  the  abscissa,  of  the  point  pi ,  and  that  x%  the  abscissa  of  the  point 
p 2.  (Nearly,  *2  <s  a  function  of  t.  The  corresponding  ordinate  ?/2  is  obtained  from  the 
line  equation  as  y 2  =  tx 2  +  /  and  is 

2t 


Note  that  is  the  familiar  parameterization  of  the  unit  circle. 

The  method  just  illustrated  can  be  thought  of  as  the  following  procedure.  Let  f(x,y)  =  0 
be  the  curve  to  lie  parameterized. 

Simple  Parameterization  Algorithm  1 

1.  Pick  a  point  p%  on  the  curve  /,  and  consider  all  lines  through  p\,  indexed 
by  the  parameter  t.  These  lines  form  a  pencil. 

2.  For  each  line  of  the  pencil,  compute  the  “other”  intersection  point  with 
/,  expressing  its  coordinates  as  functions  of  t.  These  functions  are  the 
parameterization. 

This  simple  method  works  unchanged  for  every  conic  section,  but  not  necessarily  for  curves 
of  higher  algebraic  degree.  The  reason  it  works  for  conics  is  a  consequence  of  Bezout’s 
Theorem  that  states  that  a  line  intersects  a  conic  in  just  two  points.  With  one  of  these  two 
points  fixed  as  pi,  there  is  just  one  additional  point,  and  that  point  is  uniquely  associated 
with  t.  If  /  were  a  cubic  curve,  a  line  would  intersect  /  in  general  in  three  points.  Fixing 

2  As  a  fine  point,  observe  that  the  line  *+1=0  also  contains  the  point  pi  =  (-1,0)  and  intersects  the 
circle  only  in  jm.  Here,  pi  counts  as  a  double  point,  since  the  line  is  tangent  to  the  circle.  The  point  pi  is 
therefore  the  additional  point  in  which  this  line  intersects  the  circle,  and  corresponds  to  t  =  ±00.  Only  a 
projective  parameterization  reaches  this  point  from  a  finite  parameter  value. 


f) 


one  of  them  on  the  curve,  would  leave  two  additional  intersection  points.  p2  and  P3,  both 
associated  with  the  same  value  of  t.  and  so  the  curve  would  not  be  properly  parameterized 
by  our  algorithm. 

To  make  the  method  work  for  cubics  requires  choosing  a  special  point:  that  is,  a  point 
such  that  lines  through  it  intersect  the  cubic  in  just  one  additional  point.  Such  a  special 
point  exists  only  on  singular  cubics.  and  is  actually  the  singular  point  itself.  Nonsingular 
cubics  do  not  have  such  a  point,  and  cannot  be  so  parameterized.  Indeed,  it  can  be  proved 
that  a  nonsingular  cubic  cannot  be  parameterized  at  all. 

Modified  in  this  way,  choosing  a  special  curve  point  pi .  will  our  simple  method  of  param¬ 
eterizing  curves  work  for  all  those  curves  that  possess  a  parametric  form?  Not  unless  other 
modifications  are  made  that  we  do  not  discuss.  But  there  is  an  interesting  class  of  curves 
that  can  be  so  parameterized,  namely  the  class  of  monoids,  to  be  discussed  later. 

Parameterization  of  Conics 

Conic  sections  are  parameterized  by  the,  method  described  before,  using  any  curve  point. 
However,  it  may  be  inconvenient  to  pick  an  arbitrary  point  on  the  conic.  The  computations 
simplify  if  the  point  is  picked  such  that  it  is  the  origin,  or  else  at  infinity,  in  a  principal 
direction.3  This  may  require  a  change  of  coordinates. 

Let  the  conic  equation  be 

on®2  +  2au.ry  +  a22'/2  +  2a  13®  +  2a23</  +  «33  =  0 
We  first  look  for  a  point  “at  infinity’’  by  computing  the  roots  of  the  quadratic  form 

on®2  +  fl22.'/2  +  2oi2.fi/  =  0 


The  roots  are  given  by 

X  =  —012  ±  -  «II«22  y  =  on 

They  are  either  two  real  roots,  possibly  equal,  or  they  are  conjugate  complex.  In  the  complex 
case,  the  conic  does  not  have  points  at  infinity,  and  is  an  ellipse  or  a  circle.  In  that  case  we 
must  find  a  point  at  finite  distance,  and  use  thereafter  the  method  for  monoids  described 
later.  Otherwise  we  have  a  point  at  infinity,  and  proceed  as  follows. 

Let  (k,i>)  be  one  of  the  two  real  roots.  W'e  substitute 

x  =  ®i  +  «(/i  y  =  vyi 

in  the  conic  equation.  It  is  easy  to  see  that  the  resulting  equation  is  of  the  form 

y\(cxi  +  d)  +  q(x  1)  =  0 

JSee  [1],  or  [13],  p.  177, 
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where  q(x\)  is  a  quadratic  polynomial  in  :rt  alone.  This  conic  is  parameterized  by 


*i(t)  =  /  !h  (t)  = 

ct  +  d 

Because  of  the  substitution,  therefore,  the  original  conic  is  parameterized  by 

x(t)  =  t.  +  uyi(t)  y(t.)  =  vtjiit) 

There  is  another  method  for  parameterizing  conics,  based  on  linear  algebra,  due  to  Jacobi 
in  the  previous  century.  We  sketch  it  briefly:  see  [13]  p.  170  for  details.  The  advantage  of 
Jacobis  algorithm  is  that  it  generalizes  directly  to  quadratic  surfaces,  and  that  it  is  well 
behaved  numerically. 

The  conic  equation. 


«U&'2  +  2*y  +  ft22.1T  +  2«i3.i:  +  2^23 y  +  «33  —  0 


can  be  written  as  the  bilinear  form 


(*  y  i) 


«n  a  12  a j3 

ft  1 2  ft  22  ft  23 

\  « 13  «23  ft.33  / 


W ,  \ 


y 

V1/ 


=  o 


An  iterative  algorithm  can  be  devised  that  applies  rotations  to  the  coefficient  matrix  and 
diagonalizes  it.  After  diagonalization.  t  he  resulting  conic  has  a  standard  parameterization, 
of  the  form'1 


x(t) 


/n 


i  -  t2 

i  +  /2 


y(t) 


/‘2 


2 1 

l  + t2 


The  numbers  m  and  /«2  depend  on  the  entries  bkk  of  the  diagonalized  matrix;  that  is, 
/q  =  \/IW  and  /t2  =  \/| ^22 / ^33 1 •  By  applying  the  inverse  rotations,  this  standard 
parameterization  is  mapped  to  a  parameterization  of  the  original  conic. 


Parameterization  of  Cubic  Curves 

The  difficulty  of  parameterizing  a  cubic  curve  is  to  find  out  whether  the  curve  has  a  singular 
point,  and  if  so,  where  it  is.  The  singularity  could  be  apparent  and  at  the  origin,  as  in  the 
cuspidal  cubic 

y2  -  P  =  0 

4  More  precisely,  the  parameterization  depends  on  the  signs  of  the  diagonal  entries.  We  assume  here  that 
(he  first  two  diagonal  entries  have  equal  sign  opposite  to  the  sign  of  the  third  entry. 


that  is  parameterized  by 


x(t)  =  t 2  y(t)  =  /3 

The  singularity  could  be  at  infinity,  as  in  the  case  of 

y  -  a:3  +  .r2  -  .1:  =  0 

parameterized  by 

x  =  /  y  =  /3  - 11  +  < 

But  the  cubic  could  also  be  nonsingular  and  have  no  parameterization,  as  in  the  case  of 

if  -  .r3  +  .r  =  0 

The  following  algorithm  from  [2]  solves  the  parameterization  problem  of  cubic  curves.  De¬ 
spite  the  relative  complexity,  the  algorithm  is  related  to  the  simple  parameterization  proce¬ 
dure  discussed  in  geometric  terms  before.  We  illustrate  the  method  with  an  example  from 
[i:i]  p.  181. 

We  are  given  the  cubic 

/  =  28y3  +  2  bxf  +  7  a \ly  +  f/2  +  28  y2  +  Kixy  +  7  y  +  'Ax/ 2 

Its  degree  form  consists  of  all  cubic  terms.  The  first  step  is  to  find  a  real  root  of  the 
degree  form,  so  as  to  apply  a  coordinate  transformation  to  /  that  eliminates  the  y3-term. 
.Since  the  degree  form  is  cubic,  it  has  at  least  one  real  root.  Here,  (-2,1)  is  a  root  of 
2Sf  +  26.t  if  +  7  x*y  +  :k3/2. 

In  the  second  step,  we  use  the  root  (u.v)  to  perform  the  substitution 

x  =  .Cl  +  uiji  y  =  vy{ 

Note  t hat  this  substitution  is  exactly  as  in  the  conic  parameterization.  In  our  example,  we 
obtain 

4(*i  -  1  )y?  +  4(x{  +  4.i!  1  +  l)yi  +  (.t?  +  3.x, )/2  =  0 
The  transformed  curve  equation  has  the  form 

M®i)0i  +  h(x\)yi  +  =  0 

where  li\  is  a  linear  polynomial,  I12  a.  quadratic  polynomial,  and  /J3  a  cubic  polynomial.  In 
step  3,  we  multiply  with  h\  to  bring  the  equation  into  the  form 

[h\($\)yi  +  1)]2  +  =  0 
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where  h4  =  h\h^  -  Note  that  /i,(  has  degree  4.  So.  we  have 

4(.ti  -  l)Vi  +  4(®i  -  L)(.x’i  4-  4.ri  +  I ).yi  4-  (.r|  -  x?  +  3®‘f  -  3®i)/2  =  0 
which  can  be  rewritten  as 

[2(xi  —  L)(/i  4-  ( x{  4-  4xi  +  1  )]2  -  (x'i  +  4xj  4-  l)2  +  (xj  —  x'j  4-  3xf  —  3xj)/2  =  0 

In  Step  4.  we  substitute  y\  for  the  quadratic  form  [h\(x\)y\  4-  /i2(®i))2-  In  the  example, 
i/2  =  2(*i  -  I ) J/t  4-  ( xf  +  4xi  +  1),  whence 

y\  —  ( x  *  +  1 7  Xj  4-  33xf  +  19xi  4-  2)/2 

The  righthand  side  of  this  equation  is  a  polynomial  of  degree  up  to  4.  It  can  be  proved 
that  the  original  cubic  is  nonsingular  if  and  only  if  the  righthand  side  is  of  degree  3  or  4  and 
has  no  multiple  roots.  Note  that  a  nonsingular  cubic  cannot  be  parameterized.  Otherwise, 
the  multiple  root  A  can  be  used  to  transform  the  equation  into  a  quadratic  equation,  by  the 
substitution 


X[  -  A 

In  Step  5.  therefore,  we  investigate  the  roots  of  the  righthand  side.  If  a  double  root  is  present, 
we  make  this  substitution.  In  the  example,  the  righthand  side  has  the  double  root  x\  =  -I, 
so  we  set  f/3  =  yi/(x\  +  1)  and  obtain 

2’A  —  ;,:i  T  15® i  4-  2 

If  Step  5  succeeds,  the  substitution  transforms  the  equation  into  a  quadratic  one,  and 
this  quadratic  equation  is  parameterized  in  Step  6  using  the  methods  discussed  before. 
Thereafter,  the  various  substitutions  are  inverted,  transforming  the  conic  parameterization 
to  a  parameterization  of  the  original  cubic.  In  our  example,  the  conic  is  parameterized  by 

/JV  t2-  2  #JI  /3+15<  +  2 

•Cl (t)  =  -  03  0  =  — 7= - 

2 t  +  15  v2(2t  +  15) 

Since  1/3  =  .1/2 / ( t  +  1),  we  thus  obtain 

(t2  4-  15/  4-  2)(/2  4-  2/  +  13) 

2/2  v/2(  2/4-15)2 

Now  1/2  =  2(®i  -  l)y\  +  ®2  4-  4aq  4-  1,  from  which  we  obtain 

_  (V2+l)fi  +  (8y/2+  17)/3  4-  (60y/2  4-  45)/'2  4-  (44y/2  4-  199 )<  4-  ( 109y/2  4-  26) 
yi  “  v/2(4t:}  4-  22 12  -  128*  -  510) 

Note  that  we  can  cancel  (t  -  14-  3\/2),  whence 

_  (y/2  4-  l)/34-(f)y/2  4-  12)/2  4-  (30\/2  4-  21)*  4-  ( ll\/2  4-  40) 

>Jl  ~  >/2{ -  ( 12v/2  -  26)/  -  (90\/2  4-  30)) 

From  this  parameterization,  we  finally  obtain  the  parameterization  of  the  original  cubic  /. 


y  \  =  - 
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Parameterizing  Monoids 

A  curve  /  of  degree  n  with  a  point  of  multiplicity  n  —  1  is  a  monoid.  Every  conic  section  is 
a  monoid,  and  every  singular  cubic  curve  is  a  monoid.  An  analogous  definition  for  monoidal 
surfaces  is  discussed  later. 

Monoids  are  especially  easy  to  parameterize,  provided  we  know  where  the  (n  -  l)-fold 
point  is  and  have  brought  it  to  the  origin.  (17):  or  else  if  the  point  is  at  infinity,  in  a  principal 
direction.  We  explain  the  parameterization  method  assuming  the  singularity  is  at  the  origin. 
Examples  of  monoids  in  this  form  include  the  circle 

.i;2  +  y2  -  2.c  =  0 

through  the  origin,  the  cuspidal  cubic. 

•t'3  -  if  =  0 

and  the  alpha  curve. 

.r3  +  .r2  -  y 2  =  0 

both  with  the  singularity  at  the  origin. 

When  the  [n  -  l)-fold  point  is  at  the  origin,  the  implicit  monoid  equation  is 

hn(x.y)  -  h„-i(x,y)  =  0 

where  /?„  has  only  terms  of  degree  n.  and  /?„_  i  has  only  terms  of  degree  n- 1.  This  is  readily 
verified  in  the  three  examples. 

It  is  easy  to  see  that  the  parameterization  of  the  monoid  is  given  by 


/»n— I  (•■>••  t) 

Vj  ■  . .  ■—  .i 


id-M) 


£  hn-l(s* t) 

hn(s,t) 


This  is  a  projective  parameterization  that  is  changed  to  the  normal  parametric  form  by 
either  setting  s  =  1  or  setting  t  =  1. 

The  monoid  parameterization  is  derived  by  considering  a  pencil  of  lines  through  the 
origin.  Every  line  can  be  expressed  parametrically  as 


:e(  A)  =  .sA  ij(  A)  =  tX 

and  is  determined  by  a  unique  ratio  s  :  t.  With  s  =  1,  we  obtain  the  usual  form 

H  =  tx 

Now  each  line  intersects  the  monoid  in  one  additional  point,  and  this  point  is  therefore 
uniquely  associated  with  the  ratio  .«  :  t  of  the  line. 


II 


In  the  example  of  the  alpha  curve,  :r3  +  .r2  -  y2  =  0.  we  have 

hi(x.y)  =  P 

h{x,y)  =  r-.r 

Its  (projective)  parameterization  is  therefore 

t2  -  s2 

x(s,t)  =  $ - - - 

With  s  =  I  we  obtain  the  usual  parameterization 

x(t.)  =  l1-  1  .(/(/.)  =  /3  - 12 

Because  the  parameterization  is  so  easy  to  find,  monoids  have  also  been  called  dual  forms 
in  the  CAGD  literature;  e.g.,  [17]. 

2.2  Algebraic  Surfaces 

The  parameterization  of  implicit  algebraic  surfaces  is  much  more  complicated  than  curve 
parameterization.  For  one,  a  characterization  of  when  an  implicit  surface  has  a  parametric 
form  is  technically  quite  complicated,  and  is  not  readily  explained  in  geometric  terms.  There 
is  no  known  general  algorithm  for  determining  whether  a  given  implicit  surface  can  be 
parameterized,  and  if  so,  how.  Fortunately,  monoidal  surfaces  can  be  parameterized  in  a 
very  simple  manner  with  a  clear  geometric  intuition,  and  these  surfaces  include  all  quadrics. 

Monoidal  Surfaces 

An  algebraic  surface  f(x,y ,r)  =  0  of  degree  »  that  has  an  (n  -  l)-fold  point  is  a  monoidal 
surface.  Monoidal  surfaces  include  all  quadrics,  cubic  surfaces  with  a  double  point,  and 
Steiner  surfaces,  [17]. 

Bezout's  theorem  states  that  a  line  intersects  a  surface  of  degree  n  in  exactly  n  points, 
assuming  we  admit  complex  intersections  and  intersections  at  infinity,  and  account  for  in¬ 
tersection  multiplicity.  It  follows  that  a  monoid  could  be  parameterized  by  an  extension  of 
the  simple  parameterization  algorithm  described  before: 

Simple  Parameterization  Algorithm  2 

1.  Let  pi  be  a  point  on  the  monoidal  surface  of  multiplicity  n  - 1.  and  consider 
all  lines  through  pi.  Each  line  is  determined  by  a  pair  (s,t)  of  slopes  in  two 
principal  directions,  or.  alternatively,  by  a  unique  ratio  r  :  s  :  t  of  direction 
cosines. 


y(s.l)  =  t 


1 2-s2 
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2.  Determine  (x($,t),  y(s  ,/),  z(s.t )),  the  additional  intersection  point  of  each 
line  with  the  monoidal  surface,  as  function  of  s  and  t,  thereby  deriving  a 
parameterization. 

Algebraically,  things  are  extremely  simple  when  the  singular  point  is  at  the  origin,  for 
then  the  implicit  equation  has  the  form 

hn(x.y,z)  -  hn-i(x,y,z)  =  0 

where  hn  has  only  terms  of  degree  n.  and  i  has  only  terms  of  degree  n  -  1.  The 
parameterization  is  then 

,  ,, 

x(r.s.t)  =  v - 

hn(r.s.t) 

y(r,s.t)  =  *— - 

hn{r.s,t) 

,  „  ,  />„_i(r..s./) 

z(i\sJ.)  =  I- - 

hn(r,$,t) 

This  is  a  projective  parameterization  that  is  changed  to  the  normal  parametric  form  by 
setting  one  of  the  parameters  r,  s,  or  t  to  1. 

As  an  example,  consider  the  sphere  with  radius  1  and  center  (1,0,0) 

x2  +  u2  +  z2  -  2x  =  0 


y(r,s.t)  =  s 


z(i\s.t.)  =  I 


It  is  parameterized  by 


Setting  r  =  1,  we  obtain 


x{r,s,t)  =  - - - - - 

r2  +  s2  +  t2 

2  rs 

y(r,s.L)  =  - - - - 

r2  +  s2  +  t2 

=  2rt — - 

r2  +  s2  +  t2 


x(s,t)  = 


y(s,t)  = 


Z(Syt)  — 


1  +  52  +  t2 

2  a 

1  +  s2  +  t2 
2 1 

i  +  s2  +  t2 
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Quadric  Surfaces 

Quadrics  can  be  parameterized  as  monoidal  surfaces.  Any  point  on  the  quadric  suffices.  The 
algebraic  method  discussed  for  conic  sections  can  be  generalized  [i]:  First,  we  find  a  point 
on  the  quadric  at  infinity,  next,  we  change  coordinates  to  move  this  point  into  a  special 
position,  and  then  we  pick  up  a  standard  parameterization  that  is  structurally  like  the  one 
given  before,  for  monoidal  surfaces  whose  singularity  is  at  the  origin.  The  inverse  coordinate 
transformation,  finally,  maps  the  standard  parameterization  to  a  parameterization  of  the 
original  quadric  surface. 

Jacobi’s  algorithm  can  also  be  used  for  quadric  surfaces.  Here,  the  coefficients  of  the 
quadric  surface 


«u.r  +  2<i\ixy  +  ‘la^xz  +  2aux  +  «22  y2  +  ’2«23.1/~  +  'la^y  +  n^z2  +  2a^z  +  «,n  =  0 
are  written  as  the  symmetric  matrix 


'  "11 

«I2 

"13 

"It  | 

a  12 

"23 

"21 

"13 

"23 

"33 

"34 

\  "14 

"24 

"34 

"44  J 

and  this  matrix  is  diagonalized  using  the  usual  rotations.  From  the  diagonal  form,  a  standard 
parameterization  is  determined  that  is  mapped  back  to  a  parameterization  of  the  original 
quadric  by  the  inverse  rotations.  The  standard  parameterization  depends  on  the  signs  of  the 
diagonal  entries  and  on  the  rank  of  the  matrix.  For  details  see  [13],  p.  180. 

3  Implicitizing  Rational  Curves  and  Surfaces 

Existence  of  the  Implicit  Equation 

Given  a  parametric  curve  with  rational  coordinate  functions  (/<i(s).A2(s)),  its  implicit  equa¬ 
tion  should  be  a  polynomial  f(x,y)  -  0  such  that  f(h\[s)<h2(s))  =  0.  Moreover,  the  degree 
of  /  should  be  as  small  as  possible.  To  understand  the  existence  of  the  implicit  equation,  it 
is  useful  to  learn  about  transcendental  field  extensions. 

Let  K  be  a  subfield  of  a  larger  field  E?  Consider  a  set  S  of  elements  in  E  that  are  not  in 
K.  The  extension  field  K'  =  K(S)  generated  by  S  is  the  smallest  subfield  of  E  containing 
S  and  A'.  If  every  element  of  5  is  the  root  of  a  polynomial  f(x)  with  coefficients  in  K,  then 
K'  is  an  algebraic  extension.  Otherwise.  K'  is  a  transcendental  extension.  For  example,  let 
A'  be  the  field  of  rational  numbers,  E  the  field  of  real  numbers,  and  let  S  contain  only  the 

'A  field  is  a  set  with  addition,  subtraction,  multiplication,  and  division.  For  the  properties  that  must  be 
satisfied  see  [19]. 
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number  7 r.  Then  the  extension  field  K(ir)  is  a  transcendental  extension.  On  the  other  hand, 
with  S  =  {\/2}  we  obtain  the  algebraic  extension  field  A‘(\/2).  since  \/2  is  a  root  of  x 2  -  2. 
See  also  [13],  Chapter  7. 

Consider  a  transcendental  extension  of  K  by  a  finite  set  of  elements  of  E,  say  S  = 

. t-'r } ■  It  is  then  known  that  every  element  in  the  extension  field  can  be  expressed 

as  a  rational  polynomial  expression  in  the  a:,-.  This  permits  us  to  consider  purely  symbolic 
extensions  of  K  in  which  the  elements  of  S  are  symbols.  The  resulting  extension  field  is  the 

rational  function  field  of  A*,  and  is  denoted  K(x\ . r„ ) . 

The  cardinality  of  S  is  called  the  transcendency  degree  of  the  extension,  and  a  theorem 
from  algebra  asserts  that  the  transcendency  degree  is  unique.  This  implies  that  if  uj,...,  ur+ 1 
are  arbitrary  elements  in  the  extension  field,  then  they  must  satisfy  a  polynomial  equation 
«r+i)  =  0  where  the  coefficients  of  /  are  in  A*. 

We  explain  why  the  implicit  equation  of  a  rational  curve  exists,  following  [9].  Consider 
the  rational  function  field  C(s)  of  the  complex  numbers  C  of  the  transcendental  element  s. 
Because  the  curve  is  given  by 

x  =  h\(s) 

!J  = 

we  should  think  of  x  and  y  as  elements  of  C(s).  But  the  transcendency  degree  of  C(s)  is  1, 
so  x  and  y  must  satisfy  a  polynomial  equation  f(x,y)  =  0.  The  simplest  such  equation  is 
the  implicit  equation  of  the  rational  curve. 

The  same  argument  shows  the  existence  of  an  implicit  equation  for  a  rational  surface: 
Let 

x  =  lifis.t) 
y  =  hfis,  t) 

:  =  li^(sU) 

be  the  parametric  surface.  Then  x.  y,  and  c  are  elements  of  the  field  C(s,t),  of  transcendency 
degree  2.  So.  they  must  satisfy  a  polynomial  equation  f(x,y,z)  =  0. 

Implicit  Forms  Vs.  Base  Points 

We  usually  think  of  a  parametric  curve  as  a  mapping  from  the  line  to  a  plane  curve,  and 
of  a  parametric  surface  as  a  mapping  from  the  plane  to  a  surface  in  3-space.  Suppose  we 
have  started  with  a  parametric  surface,  found  its  implicit  equation,  and  compare  the  point 
sets  defined  by  each.  Then  it  turns  out  that  every  point  on  the  parametric  surface  also  lies 
on  the  implicit  surface,  but  not  necessarily  vice-versa.  Consider  as  example  the  parametric 
surface 

(st.st\s2) 
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The  implicit  equation  of  the  surface  is 


Since 

(*/)4-(*i2)V)sO 

every  point  of  the  parametric  surface  is  also  on  the  implicit  surface.  Now  the  points 

(0. u. 0) 

are  on  the  implicit  surface.  But  when  u  ^  0.  those  points  are  not  on  the  parametric  surface, 
because  2  =  0  forces  s  =  0,  and  therefore  also  y  =  0.  It  follows  that  the  parametric  surface 
could  be  a  proper  subset  of  the  implicit  surface. 

In  general,  it  is  known  that  the  implicit  surface  may  contain  finitely  many  isolated  points 
and  curves  that  are  ■‘missing"  in  the  parametric  form  of  the  surface.  Some  of  these  points 
can  be  recovered  by  changing  the  parameterization  to  a  projective  one.  We  did  this  for  the 
unit  circle,  in  order  to  reach  the  point  (-1.0),  and  for  the  unit  sphere,  to  reach  (-1,0,0). 
But  not  all  points  can  be  so  recovered,  and.  in  particular,  not  the  points  missing  in  the 
example  just  shown. 

Consider  the  rational  surface 

Ii2{s,t)  h9(«,t)\ 

\ho(s,t)'  h0(s,t)'  ho(s.t)) 

The  polynomials  hy  define  four  curves  on  the  (^,<)-plane.  A  base  point  is  a  common  in¬ 
tersection  (iio.fo)  of  the  four  curves.  Such  parameter  values  (.?o,<o)  Jo  not  define  a  surface 
point.  Chiouh  [9]  discusses  the  relationship  between  base  points  and  “missing  points”  on 
the  parametric  surface. 

Sylvester’s  Resultant 

Sylvester’s  resultant  is  a  simple  method  for  eliminating  a  variable  from  two  algebraic  equa¬ 
tions,  and  it  can  be  used  to  find  the  implicit  equation  of  a  parametric  curve  or  surface.  Given 
two  polynomials 

f(x)  =  «„in  +  an_|.cn-1  +  •••  +  «o 
g(x)  -  bmxm  +  bm-i.  cm_1  +  •  •  •  +  &o 


it  can  be  shown  that  /  and  g  have  a  common  root  iff  the  (m  +  n)  x  (m  +  n)  determinant 


(In 

rt„_l 

flo 

0 

...  o 

0 

«n 

«i 

«o 

...  0 

0 

0 

(In 

«n- 1 

•  •  •  do 

bm 

bm — 1 

bo 

0 

...  o 

0 

bm 

bi 

bo 

...  0 

0 

.  .  . 

0 

bm 

bm- 1 

...  bo 

is  zero:  see  [19],  section  130.  The  determinant  is  the  Sylvester  resultant.  Here,  the  a*  and 
bk  are  assumed  to  be  numbers. 

The  Sylvester  resultant  can  be  used  for  two  multivariate  equations:  that  is,  under  the 
assumption  that  the  a *  and  b k  are  polynomials  in  the  remaining  variables.  In  this  case, 
the  Sylvester  resultant  R  is  a  polynomial  in  these  remaining  variables.  The  solutions  of  the 
equation  R  =  0  either  extend  to  a  common  solution  of  the  original  equations  /  =  0,  g  =  0,  or 
are  a  solution  of  the  system  an  =  0.  bm  -  0.  or  else  are  a  common  solution  of  the  coefficient 
polynomials  of  /  or  of  g. 

The  Sylvester  resultant  can  be  used  for  implicitizing  parametric  curves  and  surfaces, 
subject  to  certain  limitations.  Briefly,  the  rational  curve 

(hi{s) 

V»o  (*)’/»o(*)/ 

is  considered  as  the  intersection  of  two  surfaces  in  3-space  given  by 

xh0(s)  -  li\(s)  =  0 

yho(s)  -  h($)  =  0 

Elimination  of  s  using  the  Sylvester  resultant  delivers  the  implicit  equation,  or  a  multiple  of 
the  implicit  equation  due  to  the  possibility  that  the  lead  coefficients  vanish. 

Parametric  surfaces  could  be  implicitized  similarly,  forming  the  equations 

xho(s,t)  -  h\(sj)  =  0 
yho(sj)  -  /j2(s,/)  =  0 
zli0(s,t)  -  hz{$J)  =  0 

and  eliminating  first  s,  say,  and  then  t.  Alternatively,  a  number  of  other  resultant  formula¬ 
tions  have  been  proposed  for  eliminating  both  variables  at  once.  See  [9]  for  details. 


Grobner  Bases 

For  integral  rational  curves  and  surfaces.  Grobner  bases  algorithms  provide  an  alternative  for 
implicitizing  parametrics  without  the  introduction  of  extraneous  factors.  These  algorithms 
are  very  sophisticated  and  are  very  general:  see  [6.  7.  13.  15,  16].  Certain  specializations 
exist  that  improve  their  performance  significantly. 

We  sketch  some  of  the  ideas  that  go  into  Grobner  bases  algorithms.  For  a  more  detailed 
explanation  and  many  examples  and  applications  see  [13],  Chapter  7.  Consider  the  set 

. av]  of  all  multivariate  polynomials  in  x\ . vr  and  with  coefficients  in  the  field 

A*.  For  our  purposes,  we  consider  the  field  of  complex  numbers.  An  ideal  I  is  a  subset  of 

polynomials  in  A’(:i.*i, . tv]  that  is  closed  under  addition  and  subtraction,  and  also  under 

multiplication  with  other  polynomials  that  are  not  necessarily  in  I.  A  basic  theorem  asserts 
that-  we  can  find  a  finite  set  /j, ..../„  of  polynomials  in  /  such  that  every  other  polynomial 
U  of  /  can  be  written  as  an  algebraic  combination  of  the  .4;  that  is. 

(J  =  «i/i  +  (i2/z  +  •  •  •  +  «,./»*  €  A’(.ti . tv] 

There  are  many  ideal  bases,  not  necessarily  of  the  same  cardinality.  A  Grobner  basis  is  an 
ideal  basis  with  special  properties  that  permit  answering  basic  questions  about  the  ideal 
using  simple  algorithms.  Every  ideal  has  a  Grobner  basis,  and  this  basis  depends  on  certain 
orderings  of  terms.  For  example,  in  the  elimination  ordering 15  we  first  arrange  the  variables 
in  a  fixed  sequence,  say 

.l*t  -<  .1*2  •*!*••  ■*>  .tV 

and  declare  that  a  term  u  comes  earlier  in  the  ordering  than  another  term  v  provided  that 
r  contains  a  variable  that  is  later  in  the  variable  sequence  than  every  variable  occurring 
in  «.  So.  with  x  -<  y  -<  z ,  the  term  u  =  .r'V2  would  precede  the  term  v  =  xyz,  If  the 
highest  occurring  variables  in  the  two  terms  are  the  same,  then  the  degree  of  that  variable 
determines  the  order,  and  ties  are  broken  by  recursively  considering  subterms  derived  by 
deleting  the  highest  variable  from  both  terms.  Thus  a*10//2  -<  XUZ  and  x2yz2  <  xy2z 2. 

Given  a  set  of  algebraic  equations,  the  Grobner  basis  of  the  ideal  generated  by  the  occur¬ 
ring  polynomials,  with  respect  to  the  elimination  ordering,  defines  an  equivalent  system  that 
is  in  triangular  form  and  can  be  solved  much  more  easily.  The  basis  will  contain  the  implicit 
form  of  a  parametric  curve  or  surface,  provided  the  surface  is  integrally  parameterized. 

The  term  ordering  influences  the  time  required  to  construct  a  Grobner  basis.  The  elimi¬ 
nation  ordering  just  discussed  produces  a  basis  best-suited  to  many  CAGD  applications,  but 
requires  more  time  than  basis  construction  with  respect  to  certain  other  orderings.  Basis 
conversion  algorithms  exist  that  allow  first  computing  a  Grobner  basis  F  with  respect  to 
any  ordering,  and  then  post-processing  F  to  reveal  some  of  the  information  explicit  in  the 
basis  F'  with  respect  to  the  elimination  ordering.  Combined,  the  two  steps  often  are  much 
more  efficient  than  the  outright  construction  of  the  basis  F'.  The  approach  is  especially 

6also  called  lexicographic  ordering 
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appropriate  for  implicitizing  parametric  curves  and  surfaces.  For  details  see  [13],  Section 
7.8. 

Experience  with  Elimination  Algorithms 

Many  resultant  formulations  have  not  been  implemented  so  that  no  experimental  data  can 
be  cited  in  support  of  their  possible  practicality.  We  have  experimented  with  implicitization 
of  curves  and  surfaces  using 

1.  Sylvester’s  resultant, 

2.  Grobner  bases  with  the  elimination  order,  and 

3.  Basis  conversion. 

Three  surface  implicitization  problems  were  solved,  using  integral  parametric  surfaces  of 
degree  two,  three,  and  bicubic.  The  parametric  quadric  is 

x  —  3/2  4-  4-s2  .?/  —  2s  —  5/  T  4 

ij  —  (>s2  —  st  +  8/  +  < 
c  =  9s/.  +  12s  -  15/  +  34 

The  parametric  cubic  is 

x  —  —  /3  4-  3s/  +  s3  +  s 
;/  =  ts'2  -  3/  r  1 
c  =  2/3-5  s/  +  /-s3 

'Pile  bicubic  surface  is 

a:  =  3/(/  -  l)2  +  (s-  l)3 +  3s 
y  =  3s(s  -  l)2  +  /3  +  3/ 
r  =  -3s( s2  -  5s  +  5 )/3  -  3( s3  +  (is2  -  9s  +  i  )/2 
+Z(Gs3  +  9s2  -  1 8s  +  3)  -  3s(s  -  1) 

The  running  times  are  shown  in  the  table.  All  computations  were  done  on  a  Symbolics 
3050  Lisp  machine  with  16MB  main  memory  and  120  MB  virtual  memory.  Note  that  the 
hardware  speed  of  the  machine  is  less  than  one  MIP.  Methods  1  and  2  are  the  standard 
implementations  of  resultants  and  Grobner  bases  offered  by  Macsyma  414.62.  Method  3 
was  written  at  Purdue.  An  entry  oo  indicates  that  the  computation  could  not  be  completed 
due  to  insufficient  virtual  memory.  The  table  shows  clearly  that  Method  3  improves  effi¬ 
ciency  significantly,  but  overall  the  times  are  much  slower  than  one  would  require  for  routine 
applications.  Future  work  is  required  to  improve  the  situation. 
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Problem 

Method  1 

Method  2 

Method  3 

quadratic 

21 

22 

6 

cubic 

105 

oo 

315 

bicubic 

00 

oc 

10s 

Table  1:  Implicitization  Times  in  Seconds 
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5  The  Slides 


I.  Basic  Facts 

I. A.  Definitions  Recalled 
I.B.  Theorems  on  Conversion 
I.C.  Projective  Parameterizations 


Plane  Parametric  Curve 

x  =  hi{3)  y  =  h2(s) 
Parametric  Surface 

x  =  hi(s,t)  y  =  Ii2(s, t)  z  =  k3(s,t) 


The  hi;  are  polynomials  or  ratios  of  polynomials. 
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Parameter/Curve  Point  Correspondence 


The  parameterization  does  not  “reach”  the  point 
(—1,0)  unless  it  is  changed  to  a  projective  param¬ 
eterization,  by  homogenizing  the  functions  h 


s2  —  t2  2  st 

X~s2T~t a  y  -  -Tj-j? 

On  surfaces,  not  all  “missing  points”  may  be  so  re¬ 
covered. 


Parameter/Curve  Point  Correspondence 


Missing  Points  Example 


The  parametric  surface 

x  =  st  y  =  st2  z  =  s2 
has  the  implicit  form 

x4  -y2z  =  0 

This  surface  contains  the  line  x  =  z  =  0  that  is  not 
reached  by  the  parametric  form. 


Nonparameterizability  Example 

The  cubic  y2  —  xz  +  x  =  0  cannot  be  parameterized 
unless  square  root  functions  are  used. 

V  / 


Projective  Coordinates 

The  affine  point  (11,0:2)  corresponds  to  the  projec¬ 
tive  point  (\xq,  Aa;i,  Ao^),  where  xq  =  1  and  A  is  not 


This  sets  up  a  correspondence  between  lines  in  3- 
space  and  points  in  the  Cartesian  plane. 

Points  with  Xo  =  0  are  permitted,  but  not  the  point 
(0,0,0).  Such  points  are  called  points  at  infinity. 
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Algorithm 

1.  Fix  a  point  p  on  the  conic.  Consider  the  pencil 
of  lines  through  p.  Formulate  the  line  equa¬ 
tions. 

2.  Substitute  for  y  in  the  conic  equation,  solve  for 
x{t). 

3.  Use  the  line  equations  to  determine  y(t). 


Implicit  equation: 

xy  —  1  =  0 

Lines  are: 

1 

-♦o 

r-H 

+ 

1! 

Substitution  yields: 

x(x  +  l)t  —  x  —  1  =  0 

Interesting  root: 

x(t)  =  l/t 

Parameterization: 

x(t)  =  l/t 

y(t)  =  t 
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How  to  Find  a  Point  on  a  Conic 

Intersect  with  a  line  -  requires  solving  a  quadratic 
equation. 

Easiest  with  line  at  infinity,  but  that  may  yield  com¬ 
plex  points. 

Otherwise,  find  extrema  by  intersecting  with  a  par¬ 
tial. 
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Example 


Given: 
Degree  form: 
Root: 

Substitution: 


x2  -  y2  +  2x  +  y  +  4 
x2-y 2 
(-1.1) 


*  =  *i  —  2/1 


y  =  2/i 


Result:  +  2xx  +  4  —  yi(2xx  +  1)  =  0 

Parameterization:  xx  =  t 

2/i  =  (£2  +  2 1  +  4)/ (2<  +  1) 
t2  _  i  _  4 

Backtransformation:  *(<)  = - 

W  2t  +  1 


2/(0  = 


t2  +  2t  +  4 
2<+l 


Jacobi’s  Method 

1.  Write  the  conic  os  the  bilinear  form 

an  a12  aj3  ^ 

M 

(*  V  0 

Ol2  <*22  a23 

y 

=  0 

^  <*13  ®23  a33  ) 

l1 ) 

2.  Diagonalize  by  applying  rotations  B 

=  RART 

3.  Backtransform  a  standard  parameterization  of 
the  diagonal  form. 

Otherwise,  assume  without  loss  of  generality  that 
At  >0.  If  A2  is  also  positive,  then  the  parameteriza¬ 
tion  is 

...  1  -t‘ 

*M  =  ft-rrj 


where 


»(0  = 


Ml  =  v/|Al/A3[  M  2  =  v/ 1 A2/  A3 
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X. 
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A  Worked  Example 


28y3+26xy2+7x2y+x3/2+28y2+l6xy+7y+8x/2  =  0 
Degree  form  is 

28  y3  +  26  xy2  +  7x2y  +  x3/2 
and  has  the  real  root  (—2,1). 


Substitute 

x  =  xi-  2yi 

y  =  yi 

Yields 

4(xi  -  l)y?  +  4(xJ  +  4xi  +  l)yi  +  (xi  +  3xi)/2  =  0 


Multiply  with  (xi  —  1)  and  regroup,  obtaining 

(2(xi  -  l)yi  +  (*1  +  4*!  +  l))2 
-  (x2  +  4xi  +  l)2  +  (xf  -  xj  +  3xf  -  3ii)/2 

Substitute 

y2  =  2(xi  -  1)3/1  +  (x2  +  4xt  +  1) 

Result  is 

y\  =  (*{  +  17x3  +  33x?  +  19*i  +  2)/2 


Here,  -1  is  double  root  of  righthand  side.  Substitute 


X\  +  1 


Result  is 

2yl  =  x\  +  15xi+2 
which  is  parameterized  by 

*1  =  M  2/3  =  ~i 


f2  4-15*4-2 
\/2(2t+lS) 
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II.B.  Surface  Parameterization 
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Implicit  Equation  of  Monoids 

If  the  (n  —  l)-fold  point  is  at  the  origin,  the  implicit 
equation  of  the  monoid  is  of  the  form 

hn{x,y,z)  y,z)  =  0 

where  hn  has  terms  only  of  degree  n,  and  hn-i  has 
terms  only  of  degree  n  —  1 


Parameterizing  a  Monoid 

Monoid  hn(x,y,z)  -  hn-i(x,y,z)  =  0 

Parameterization  x(s,t)  =  -  — ^ 

hn{l,s,t) 

MM,*) 

z(s,t)  =  t^ izlihhH 

MM,*) 


xample 


(a:2  +  y2  +  z2)  -  2x  =  0 
has  the  parameterization 


End  of  Part  II 
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The  general  element  in  the  transcendental  extension 
K{si,...,$m)  has  the  form 

P(«l  5  •••»  ) 

U  =  - 

where  p  and  q  are  polynomials  with  coefficients  in  K. 


Algebraic  Dependence 

Let  K'  be  a  transcendental  extension  of  A',  obtained 
by  adjoining  finitely  many  s*. 

The  elements  in  K'  -  K  are  algebraically 

dependent  if  there  is  a  polynomial  q(xu  ...,xr)  with 
coefficients  in  K  such  that  the  u*  are  a  root  of  q. 

Otherwise  the  ujt  are  algebraically  independent. 


Transcendency  Degree 

The  transcendency  degree  of  K'  is  a  number  d  such 
that  any  d  +  1  elements  in  K'  —  K  are  algebraically 
dependent. 

Theorem 

The  transcendency  degree  of  K'  is  unique. 
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Example 

The  parametric  surface 

x  =  st  y  —  st2  z  =  s2 
has  the  implicit  form 

x4  —  y2z  =  0 

This  surface  contains  the  line  x  =  z  =  0  that  is 
not  on  the  parametric  surface,  except  for  the  point 
(0,0,0). 


IILB.  Sylvester’s  Resultant 
Variable  Elimination 


Given  two  polynomials 

f{x)  =  anxn  +  an_i.Tn“1  + - H  a0 

g(x)  =  bmxm  +  +  •  ••  + 


we  want  a  criterion  for  a  common  solution. 
Form  the  determinant  R 
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Intuition 

Given  a  system  of  linear  equations,  manipulations 
such  as  LU-de composition  derive  an  equivalent  linear 
system  that  is  easier  to  solve. 

Likewise,  given  a  system  of  nonlinear  equations,  a 
Grobner  basis  is  an  equivalent  system  of  nonlinear 
equations  that  is  easier  to  solve. 
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The  unique  representation  is  an  ideal... 

Also  true  for  curves,  surface  intersections,  and  so  on. 


Ideals 

An  ideal  I  is  a  set  of  polynomials  such  that 

1.  If  p  and  q  are  in  /,  then  so  is  p  -  q. 

2.  If  p  €  /,  and  q  is  any  polynomial,  then  pq  is 
also  in  I. 


Ideal  Bases 

All  ideals  are  finitely  generated;  that  is,  there  are 
polynomials 

A, /a,  M,i  fm 

in  /  such  that  every  other  polynomial  in  /  can  be 
written 

q  =  «l/l  +  «2A  +  •  •  •  Mm/m 
where  the  u,-  are  polynomials. 


A  Grobner  basis  is  an  ideal  basis  with  special  prop¬ 
erties. 

In  particular,  a  GB  with  respect  to  the  “elimination 
ordering”  allows  solving  a  system  of  nonlinear  equa¬ 
tions  in  an  especially  simple  manner. 
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The  parametric  cubic  is 


x  =  —  t3  +  3  $t  +  s3  +  $ 
y  =  ts2  -  3 1  -f  i 
z  =  2t3  —  5st  +  t  —  s3 


Implicitization  Times  in  Sec. 


Problem 

Method  1 

Method  2 

Method  3 

quadratic 

21 

22 

6 

cubic 

105 

oo 

315 

bicubic 

oo 

oo 

105 
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